Threshold singularities in the dynamic response of gapless integrable models 
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We develop a method of an asymptotically exact treatment of threshold singularities in dynamic 
response functions of gapless integrable models. The method utilizes the integrability to recast 
the original problem in terms of the low-energy properties of a certain deformed Hamiltonian. 
The deformed Hamiltonian is local, hence it can be analysed using the conventional field theory 
methods. We apply the technique to spinless fermions on a lattice with nearest-neighbors repulsion, 
and evaluate the exponent characterizing the threshold singularity in the dynamic structure factor. 
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One-dimensional (ID) interacting systems [liSill] oc- 
cupy a special place in quantum physics. Although in- 
teractions have much stronger effect in ID than in higher 
dimensions, it is often possible to evaluate observable 
quantities exactly. Besides forming the basis of our un- 
derstanding of strong correlations, ID models have long 
served as a test bed for various approximate methods. 

Many properties of interacting ID systems can be 
understood by considering integrable models [il]. The 
paradigmatic example is N spinless fermions on a lattice 
with L cites with periodic boundary conditions, 

L 

H =J2[~ ^In'^m+i - ^L+i^m + 2An,„n™+i . (1) 

m— 1 

Here rim — 4'm''Pm ^i^d A is the repulsion strength 
The experimentally relevant response function for the 
model ^ is the dynamic structure factor 



S{q,uj) 



2Tr 



J2\{f\nl\0)\'s{u;~ef + eo) (2) 



in the thermodynamic limit L — )■ oo taken at a con- 
stant filling factor = N/L < 1/2. In Eq. = 

Ek4A-, with Vfc = i-^/^E™e-™Vm, l/> is an 
eigenstate of Eq. ((!]) with energy e/, and |0) is the ground 
state with eo being the ground state energy. 

In any ID system, conservation laws restrict the sup- 
port of correlation functions in (w, q) plane. For example, 
S{q,uj) > only at w > cuo{q), see Fig. [TJ On general 
grounds, S{q,uj) is expected to exhibit a power-law sin- 
gularity at the threshold @ , 



S{q,uj) (x[uj~ uJo{q)] 



(3) 



Although exact eigenstates of integrable models can 
be constructed using the Bethe ansatz Q, evaluation of 
dynamic correlation functions is very difficult. A con- 
siderableprogress was achieved in understanding gapful 
models [Tj. However, it is still largely an open prob- 
lem in the gap less case, and, with few exceptions (see, 
e.g., [1, lloL ITl|). the threshold exponents /i are not 



known. For the model ([T]), the most complete results so 
far were obtained b y co mbining numerics with the alge- 
braic Bethe ansatz jl2l [l3| . The main limitation of this 
technique is very slow convergence towards the thermo- 
dynamic limit, which makes it very difficult to evaluate 
the exponent. 




FIG. 1: Support of the structure factor in (uo,q) plane. For 
q 7^ 0, 2fcF, the boundary of the support uj(i{q) lies below the 
straight dashed lines uii{q) = min{ug,«|g — 2fcF|}; here kp ~ 
■KU is the Fermi momentum. 

In the alternative Luttinger liquid approach [l|, , 
a ID system is described by two parameters, the sound 
velocity v — {cLuq / dq) , and the Luttinger parameter 
K characterizing the interaction strength. At q — *■ 2fci?, 
Luttinger liquid theory yields Eq. ([3]) with the exponent 



ML 



1 



(4) 



At q —f 0, however, one finds S oc 6{lu — vq). The 
discrepancy with Eq. ([3]) is due to the omission in the 
fixed-point Luttinger liquid Hamiltoni an jl , [l^ of the 
irrelevant in the RG sense operators IJ, [iJ] representing 
the spectrum nonlinearity. Indeed, for small q, most of 
the spectral weight of S{q,uj) is confined to a narrow 
interval of uj of the width Slo ^ uji — ojq about w = wi, 
while Eq. ([3]) is applicable at <C 5uj 0]. For a hnear 

spectrum 6lo — 0, which makes the regime of Eq. ([3]) 
inaccessible. For the same reason, the Luttinger liquid 
result (HI) is valid, strictly speaking, only aX w — Wq ^ 5uj, 
and the exact exponent may differ from /i^. 

In fact, the true threshold exponents often deviate from 
the Luttinger liquid theory predictions already in the low- 
est order in the interaction strength for S{q,uj) near 
q — 2k F, such deviations show up at q > 2k p. 
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In this Letter we develop a technique that allows exact 
evaluation of the exponents characterizing threshold sin- 
gularities. The technique is applicable to any correlation 
function that exhibits a threshold behavior and to any 
gapless model integrable by the Bethe ansatz. We illus- 
trate the idea of the method by working out the dynamic 
structure factor ([2]) for the model ^ as an example. 

Among various states |/) for which (/| jt-J |0) [see 
Eq. ([2|)] for a given q, one state, say, \fq), has the low- 
est energy = eg + ^oil)- Consider now a deformed 
Hamiltonian H with the following properties: 

(i) The deformed Hamiltonian H is local. 

(ii) H commutes with H, so that H and H have the 
same set of eigenstates |/). 

(iii) The states |0) and \fq) represent the doubly- 
degenerate ground state of H. 

(iv) The deformation H ^ H is continuous in the sense 
that if the state |/) has momentum q and its energy 
e/ is close to e^, then the corresponding eigenvalues 
of H satisfy ef - eq = e/ - + 0[(e/ - 6,)^]; 
similarly, e/ — eo — ^ + O [(e/ — eo)^] for states 
I/) with zero momentum. 

Once the Hamiltonian H satisfying these conditions 
has been constructed, it can be analysed using conven- 
tional field-theoretical methods [3j. In particular, condi- 
tion (i) allows one to identify the low-energy projections 
of microscopic fields with local operators in the effective 
continuum description. The coupling constants of this 
effective low-energy theory can be found by comparing 
its low-energy spectrum with that of H. (Essentially the 
same ideas are behind the Luttinger liquid description 
of the low-energy excitations of spin chains ^). Finally, 
conditions (ii)-(iv) guarantee that the structure factor 
calculated for the Hamiltonian H will have a power-law 
singularity at u —^ with the same exponent fi that char- 
acterizes the threshold behavior in the original model. 

Integrable models have an infinite number of indepen- 
dent local operators /„ commuting with H (integrals of 
motion). Thus, any Hamiltonian of the form 

H = Y, Cnin (5) 

will satisfy conditions (i) and (ii). 

For free fermions \fq) = iplpi^kp-q 1^)' ^ig. [21(a), 
and the integrals of motion have a simple form. Con- 
ditions (ii)-(iv) will be fulfilled if the single-particle 
spectrum of the deformed Hamiltonian has the shape 
sketched in Fig. [2](b). While this is not sufhcient to de- 
termine the coefficients c„ uniquely, the low-energy spec- 
trum of H is completely specified. 

In the Bethe ansatz, excitations of integrable mod- 
els such as Eq. (m are classified in terms of fermion- 
like quasiparticles and quasiholes 0. Similar to free 





FIG. 2: (a) For free fermions (A = 0) and for u < 1/2 
and q < 2kF, the state \fq) consists of a single particle- hole 
pair added to the Fermi sea: the particle is created at the 
Fermi momentum k = kp, and the hole ai k — kp ~ q. 
(b) Single-particle spectrum of the deformed Hamiltonian H: 
only states in a narrow strip of energies (shaded) about the 
Fermi level (dashed line) contribute to the structure factor 
near the threshold. 



fermions, the state \fq) corresponds to a quasiparticle 
with momentum kp and a quasihole with momentum 
kp — q added to the ground state, cf. Fig.[2](a). 

Provided that the quasiparticle spectrum of the de- 
formed Hamiltonian (see below) is similar to that shown 
in Fig. [21(b), it is obvious that the infrared fixed point 
of H corresponds to a single hole minimally coupled 
to the right- and left-moving fermions with linear spec- 
trum. We introduce bosonic fields {p± which satisfy 
\ipa{x),'^a'{y)\ = JTra sgn(a:: — y), and the field d 
which describes an infinitely heavy hole [l^ . The fixed- 
point Hamiltonian then assumes the form familiar from 
the x-ray edge singularity problem [l7j . 



H = 



'^'^Vo^^idx^af ~'^Pa{dx^a)d{x)d'^{x) . (6) 



The low-energy projection of the microscopic field V'm 
is given by ^Am = Y^^e'"'"=•''^pa{x)+ e}('"^-''')''d{x), where 
m = a; is treated as a continuous variable and the fields 
'0± are related to 1^9-1- according to 

V'a oc exp[ia(iy9Q. cosh?? — iy9_Q, sinhi?)] , e~'^^ = K. (7) 

The leading contribution to the density operator [see 
Eq. ([2])] is then given by 



(X / dx ip\,{x)d{x). 



(8) 



Evaluation of the structure factor ([2|) using Eqs. ([6|)- 
yields a power-law singularity with the exponent 



, = l-(cosh.-^|±)^(sinh. + ^^' 



(9) 



We now sketch the construction of the deformed 
Hamiltonian (O and the derivation of the coupling con- 
stants of the corresponding fixed-point Hamiltonian (|6|. 
We choose /q = J^m'^ti'^Tm so that cq in Eq. (O plays 
the role of the chemical potential. For n > 0, the inte- 
grals of motion are expressed via the derivatives of the 
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transfer matrix (trace of the monodromy matrix) t(^) 
with respect to the spectral parameter ^ |2] , 



^iTT /2—iri ' 



(10) 



The first operator in this hierarchy is proportional to the 
Hamiltonian itself: Ii — H/ sin(2?7). The next one, I2, is 
given in a closed form in [l8| . 

Consider a quasiparticle (quasihole) excitation of the 
Hamiltonian ([1]) characterized by the rapidity A [l^ . By 
construction, such excitation is an eigenstate of the de- 
formed Hamiltonian i?, see Eq. ([5]). Using properties of 
the transfer matrix [3| , one can show [20] that the corre- 
sponding eigenvalue e^(A) satisfies the equation 



(11) 

Here the Fermi rapidity Aj? is the solution to fc(Ai?) = 
kp [l3|, K{X) = dO /dX where 0{X) is the bare two-particle 
phase shift, and po(A) is the bare particle momentum; for 
the model ([T|) these are given by Q 



= i In 



sinh(A + 2iri) 



sinh(— A + 2m;) 



Po = « In 



cosh(A — irj) 



cosh(A + irj) 



(12) 

In order to satisfy the conditions (ii)-(iv) above, we 
impose additional constraints on e^(A), 



(de/dA)±A^ 



7{± Xf) = 0, 
= {de/dX)±x, 



{dI/dX)x, 

_P(A£) 

P(A,) 



0, 

{de/dX)x, 



(13) 



where Xq is the solution to k{Xq) = kp — q [l9|. The 
constraints are equivalent to five linear equations on the 
coefficients c„ in Eqs. ([5]) and (fTTjl . In order to sat- 
isfy these equations, it is sufficient to keep the first five 
integrals of motion in Eq. ([5]). 

The coupling constants of the effective fixed-point 
Hamiltonian ^ follow from the comparison of the finite- 
size spectrum of © with the exact low-energy spectrum 
of the deformed Hamiltonian H. This procedure is stan- 
dard 3 and yields 



/3± =±27ri^(± Xp\Xq), 
where F is the dressed phase shift that satisfies 01 



(14) 



dfiKiX- 



^t^)F{fi\O^^0iX^ 



<)■ (15) 



Eqs. and (fT5|) uniquely define the parameters /?±, 

and, therefore, the exponent ([9]). 

Note that (3± do not depend explicitely on the coeffi- 
cients c„. Indeed, the constraints (|13|) do not completely 
fix H but only specify its low-energy spectrum. We em- 
phasize that our construction does not rely on the model- 
specific Eq. ([T^ but is applicable to any model integrable 
by the algebraic Bethe ansatz. 



We now use Eqs. ([9]), ([M)) . and (fT5|) to evaluate the 
threshold exponent for the model ([1]). Precisely at half- 
filling Xp 00 2] and Eq. ([T5]) is solved by Fourier 
transform with the result 



F(±oo|C) = ±(k-1)/2, 



(16) 



where we used the well-known |2|] value of the Luttinger 
parameter at half-filling, k = tt/Aj] Eqs. (O, (fTi)l . 
and (|16p then yield a momentum-independent exponent 



1 



1 

K 

K 



(17) 



Comparison with Eq. Q shows that the exact exponent 
fiQ is smaller than the Luttinger liquid result /i^ (note 
that K varies between 1/2 and 1). For a weak interaction, 
/iL"Mo~Mi/2»2(A/7r)2. 




FIG. 3: Threshold exponent ^{q) for A = 0.9 and for 
f = 0.4,0.45,0.49,0.495,0.499 (bottom to top). The dashed 
horizontal lines correspond to /il and fio at half-filUng. 



Away from half-filling, Eq. ([15]) can be solved nu- 
merically. The resulting exponent is g-dependent, see 
Fig. [3l It varies from ^(0) = to n{2kp) = ^2 with 
Mo < ^J'2 < the exact value of /j,2 depends on both A 
and V [iO]. Very close to half-filling, the dependence ^{q) 
is nonmonotonic: outside narrow intervals of the width 



- 1/2 - < kp 



(18) 



near q — and q ~ 2k p, the exponent approaches a 
constant, ^{q) « /ii = const. Surprisingly, //i coincides 
with the Luttinger liquid exponent fi^ rather than with 
the exact half-filling result fiQ. 

This discrepancy originates in the peculiar behavior of 
the phase shifts near the Fermi points. Consider Eq. pS]) 
at ± Ai?| ^ i5A with SX = 1 — 2rj/Tr (this limit corre- 
sponds to q away from q = 0, 2kp). In order to find the 
phase shift -F(A|C) for A^ 3> 5X (i.e., close to half-filling) 
and Xf^Xp (i.e., close to the right Fermi point), we re- 
place 6{X — C) in the r.h.s of (jlSp by 0(oo), and extend the 
integration in the l.h.s. to —00. The resulting equation 



^(AIC)-:^ / ' diiK{X~y)F{ii\0^^e{^), (19) 
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as well as the similar equation for A « —Xf, describes 
the fractional charge function Its solution yields 



2^F(A|C) 



70, A, Af — A > 

71, A_F — A <C <5A 



(20) 



with 7o = n{K — 1) and 71 = tt{k — In other 

words, the limits A — > ± Af and Xp 00 for the phase 
shift F(A|C) do not commute. If the limit A ^ ± A^? is 
taken first, and the resulting phase shifts are substituted 
into Eq. p^ . one finds /3± = 7i- Eq- ® then yields 
fii ~ fi^. However, by taking first the limit Xp ~^ 00 
(i.e., V —^ 1/2), one would find /3± — 70 and fix — fiQ. 

The noncommutativity of limits has observable conse- 
quences. Indeed, F(A|C) characterizes the strength of the 
interaction between a quasihole at rapidity C and a quasi- 
particle at rapidity A. According to Eq. the phase 

shift at A ~ ±Ai? changes with A on the scale 6X Xp- 
In the momentum space, this corresponds to narrow in- 
tervals of the width 6k = 2ttp{Xf)SX ~ (1/2 — v) near 
k ^ ±kp (here we used the well-known Q result for p{X) 
near half- filling) . States within or outside these inter- 
vals interact with the quasihole with coupling constants 
f3± 71 or P± « 7o, respectively. 

As the filling factor approaches 1/2, the interval 6k 
collapses. For a finite-size system close to half-filling, 6k 
will eventually become compatible with the momentum 
quantum ~ 1/L. In this limit, the threshold behavior is 
dominated by states outside the interval Sk. Accordingly, 
the exponent fii that characterizes the threshold singu- 
larity away from q = 0,2kF exhibits a crossover from 
/xi = /LiL at 1 ^ (5fc 3> 1/L to fii = Ho at (5fc ^ 1/L. 

In the recent study (211 the exponent fi at half-filling 
was found to be equal to (Xl- Our consideration shows 
that fi indeed approaches this value when v 1/2. How- 
ever, because the limits v —^ 1/2 and w — > do not com- 
mute, the region of applicability of the result fi — fiL 
limited to uj — uiq v6k. Precisely at half- filling 6k — 0, 
and the exponent is given by Eq. p7p instead. 

It should be mentioned that the two-spinon contribu- 
tion to the structure factor has a square-root singularity 
at w ^ Wo [22] . This result was obtained by approaching 
A = 1 from the gapful side of the transition A > 1. We 
found /zo < 1/2 in the gapless regime A < 1, see Eq. p7)) 
above. The discrepancy suggests that the threshold ex- 
ponent Ho{^) has a discontinuity at A = 1. 

Finally, for 6k > and when q approaches either 
or 2kp, the situation is complicated by the competition 
between two small energy scales, v6k and 610 (see above). 
The behavior of the structure factor in this regime will be 
discussed in details elsewhere bd]. At small g, it agrees 
with the first-order result of [6|: /x(g) ~ {A/6k)q for 
< (5fc ^ 1 and /io ~ /^l ^ A for 6k = 0; the two 
exponents merge at q ^ 6k. 

To conclude, we proposed a method of evaluating the 



exponents characterizing threshold singularities in the 
dynamic response functions of gapless integrable mod- 
els. Application of the method to the dynamic structure 
factor of ID spinless fermions on a lattice revealed un- 
expected complexity in the dependence of the threshold 
exponent on the system parameters near half-filling. 
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